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We show that by a suitable choice of a time dependent Hamiltonian, Deutsch's algorithm can 
be implemented by an adiabatic quantum computer. We extend our analysis to the Deutsch-Jozsa 
problem and estimate the required running time for both global and local adiabatic evolutions. 



Quantum computation and quantum information theory has attracted a great deal of attention in recent times. 
Inherently quantum mechanical systems can in principle be used to implement a wide variety of computational 
algorithms wth enhanced efficiency [1-3]. The principle of superposition in quantum mechanics, according to which a 
system can be in a linearly superposed state of more than one eigenstate, is the key to this increased efficiency. One 
of the first algorithms that was first proposed in this context is Deutsch's algorithm [4]. 

In this, one would like to determine whether a function 



Ordinarily, one has to determine both /(O) and /(I) to infer the nature of the function, since the knowledge of one 
does not shed light on the value of the other. However, it was shown that by applying a certain sequence of unitary 
operators ('gates') on a given initial quantum mechanical state, and then making just one measurement on the final 
state, the nature of the function / can be determined [4] . 

Recently, a new framework of quantum computation has been proposed, in which the series of gates referred to 
above is entirely replaced by a Hamiltonian which changes continuously with time. The Hamiltonian is so chosen 
that the state of the system is its ground state at all times (although the ground state itself is time dependent), and 
the system slowly evolves to a desired final state [5]. Several applications of this have been considered [6]. Using 
this framework, it was shown that Grover's search algorithm can be efficiently implemented [7]. In this paper, we 
show that Deutsch's algorithm can be implemented as well, by choosing a suitable initial state and a Hamiltonian 
which evolves that state. Then a single measurement of the final state suffices to determine whether the function / 
is constant or balanced. Finally, we show that the results can be extended to the Deutsch-Jozsa algorithm involving 
n-qubits. 

Let us begin with a 2-level system, e.g. a spin 1/2 particle, with the basis kets {|0), |1)}. We define the 'initial' and 
'final' Hamiltonians Hq, Hi respectively as: 



/: {0,1} ^{0,1} 



is constant or balanced, i.e. whether /(O) = /(I) or /(O) ^ /(I) using a quantum computer. 
The four possible outcomes of / are: 



/(O) = /(I) (constant) 

/(O) = /(I) = 1 (constant) 

/(O) = , /(I) = 1 (balanced) 

/(O) = 1 , /(I) = (balanced) 



Ha 
Hi 



/-|^o)(^o| 
/-|V'i)(^i| 



(1) 
(2) 
(3) 



where the initial and final state vectors are 



;iven respectively by: 



^(10) + 11)) 
a\0)+m 



(4) 



IV'i) 



(5) 
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with 



a 



(3 



(_l)/(0) + (_l)/(l) 

(_l)/(o) _ 



(6) 
(7) 



Note that the Hamiltonians in the above are implicitly given in terms of some fundamental energy scale, E, say, whose 
value is determined by the physical system used to construct the states. This energy scale has a natural time scale 
associated with it, namely T = U/E, which will play an important role later. The following relations will also prove 
useful: 



a + p = 1 

a^ = a , 0^ = (3 
& a/3 = 



(8) 
(9) 
(10) 



Thus, when / is a constant, a = 1, /? = 0, and vice-versa. We assume that the system is initially in a state |V'o) and 
is evolved by the following time dependent Hamiltonian: 



H{t) = (1 - s{t)) Ho + s{t) H, 



(11) 



In general s{t) is an arbitrary function of the time t, such that s(0) = and s(l) = 1. Thus, H{0) = Hq and 
H(l) = Hi. For the present we assume that s{t) is linear in t, namely s = t/T, where T is the total time for which 
the system is evolved. 

It follows from the adiabatic theorem that at i = T, the system would be in the state I'^i), with very high probability 
(1 — e^)^ (where e is a small number), provided the evolution is adiabatic [8, 9], i.e. 



< e 



(12) 



where the time t is implicitly given in units of T. The lower bound on the evolution time T is implicitly given by the 
condition (12) as will be seen later and 



E^ 



dH \ _ ds IdH 
~dt / ~ dtXlJ 



dH 



dt 



E. 



(13) 

(14) 



E^{s) and E_{s) being the two time dependent eigenvalues of H{s), with corresponding eigenvectors \E^,s) and 
\E-,s) respectively, and 

Qmin ■■= min [E+{s) - E-{s)] 

0<s<l 

Of course, the quantities < dH/dt > and gfmin should be non- vanishing. We will show that this is indeed the case. 

Thus, under adiabatic quantum evolution, the initial wave function of the system will evolve to either |0) or |1), 
almost with certainty, and by making a single measurement of the state at the end, one can determine unambiguously 
whether the nature of the function / is. Namely, if the outcome of measurement is |0), then / is constant and if it is 
|1), then it is balanced. 

The matrix elements of H in the {|0), |1)} basis are: 



H{s) 

The corresponding eigenvalues are: 



1/2 + s(/? - 1/2) 
-1/2(1 -s) 



- 1/2(1 -s) 
1/2 + s(a- 1/2) 



1 ± Vl - 2s + 2s2 



(15) 
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Note that the eigenvalues arc independent of a and /3, and are identical to those found in [7]. This is a consequence 
of the relations (8) - (10). It follows that 

Ai;(s) = £+(«) - E_{s) = Vl - 2s + 2s2 , 
which is non-zero for all values of s and 



= AE{s = 1/2) = 



(16) 



In addition, the two orthonormal eigenvectors are: 



\E±,s) = k± 



\ [(1 - 2a)s T Vl - 2s + 2s^\ J 



where 



with which one can get: 



fc± = 2-1/2 [(1 -2s + 2s^) ± (2a - l)sVl - 2s + 2s2 



-1/2 



dt 



T 2\/l - 2s + 2s2 



(17) 



Note that this too is independent of a and /3 although \E±, s) is not. Also, note that for any value of s, the above 
quantity is non-vanishing and of order 1/T. At the final time t = T, the eigenstates are |0) and |1) respectively, and 
depending on the value of a, (3, the system evolves to one of them. 
Substituting (16) and (17) is (12), we get the following relation: 



T>1 
e 



(18) 



which gives an estimate of the time for which the initial state |V-'o) must be evolved via the adiabatic Hamiltonian 
(11) to attain an accuracy of order e of the final result. For example, if we want the final state to be the \ipi) with 
accuracy of 90%, then the minimum evolution time should be order of T w — \/0^ w 4.4, in units of T. 

A few comments arc in order here. Instead of the starting with the initial state (4), one can in general start with 
an arbitrary initial state of the form 



IV-o) =a|0) + 6|l) 



(19) 



with |ap + |6p = 1, and evolve the system with the Hamiltonian (11). The end result is expected to remain unchanged, 
since the ground state of the final Hamiltonian is still to which the system will eventually tend. If one starts 
with the state (19), then A(s) and {dH/dt) are respectively: 



A(s) = ^/l-4{a^/3 + b^a)s{l-s) 



dH_ 
dt 



ah 



T ^1 - 4(a2/J + 62q;)s(1 - s) 
Substituting (20), (21), s = 1/2 in (12) and siniphfying we get: 



T> — 



ah 



Thus depending on whether a = of a = 1 (although this value is a priori unknown), we get respectively: 



T> 



1 a 



or 



T > 



1 h 



e a 



(20) 
(21) 



(22) 



(23) 
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Thus for extreme asymmetric values of a and h (e.g. a w and 6 « 1), the evolution into the final state would either 
take place in a very short or a very long time. But as the value of a is not known, one would have wait for the greater 
of two values before making the measurement. Equivalently, for very small a or h, gmin becomes very small for some 
value of a, which is contrary to what the validity of adiabatic theorem requires. Consequently, the 'optimal' values 
for which the evolution time is independent of a, /? is given by 

a = b = ^= 
V2 

which is what we started with. 

The above procedure can be generalized to Boolean functions of the form : 

/: {0,1}" ^{0,1} 

by making use of n qubits instead of a single one [1, 10]. In accordance with the Deutsch-Jozsa problem, we assume 
that it is 'promised' that the function is either constant (i.e. all outputs arc identical) or balanced (i.e. has an equal 
number of Os and Is as outputs), and the task is to find which of the above it actually is [1, 10]. The basis states now 
are {|0), ]1), . . . , jA'' — 1)} (with N = 2"). Now we choose the normalized initial and final states to be: 



I^o) = :^EK) (24) 



i=0 



N-1 



l^i)=alO) + -^^5]lfc) 



(25) 



k=l 



with 



a = — 



xe{o,i}" 
a 



(26) 
(27) 



Once again, if f{x) is constant then a = 1 and /? = and vice- versa. Thus after the required running time, if a 
measurement of the final state yields ]0), then f{x) is constant and if it does not yield )0), then it is balanced. The 
properties (8 - 10) continue to be valid. Also, Hq,H\ and H{s) are still given by (1), (2) and (11) respectively. 
In the chosen basis, the adiabatic Hamiltonian i?(s) is now given by the following N x N matrix: 



/I 



H{s) 



l-s 
N 



— as 



l-s 
N 



l-s 
N 



N 



\ 



N 



N 



N 



_s§_ 

N-1 



N 



_s§_ 
N-1 



N 



N 



_s§_ 

N-1 

_S0_ 

N-1 



N 



N 



sfi 
N-1 

N-1 



\ N 



l-s _ sl3 
N N-1 



l-s _ s0 
N N-1 



Af N-1 / 



It can be shown that the highest eigenvalue of the above Hamiltonian is 1, which is (A'' — 2)-fold degenerate, and the 
two remaining distinct eigenvalues (both less than 1) are: 



l±^/l_l£(l_^(/3 + a(7V-l)) 



(28) 



Thus, 



!^E{s) = E+{s) - E. is) ^x/l- ^^^^ (/? + a{N - 1)) ^ , 



(29) 
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implying 



AE(s = l/2) = \/l--(/3 + a(7V-l)) 



In addition, it can be shown that: 



dH 
ds 



y/N - 1 

N 



4s(l- 



N 



■[/3 + a(iV-l)] 



Thus for s = t/T, from condition (12) it follows that (for large A'') 

T> — 



(30) 



(31) 



(32) 



which shows that the evolution time scales as N, the number of qubits. 

However, following [7] if we assume an evolution with a general s{t) we obtain the adiabaticity condition (12) that 
must be satisfied at any given instant of time t: 



ds [E+{s)-E_is) 
1(^)1 

Substituting (29) and (31), we get (here 7 = 4(/J + a{N - l))/N) : 



1 ^jin^ 

dt> - -^^77 ' 

- e N 



I [1- 



57 + 3^7]^/^ 



(33) 



(34) 



The expression on the r.h.s. can be integrated, and is equal to 



x/F^ 1 



2s- 1 



N 



y3/2 



(7-1 - l/4)Vs2-s + 7- 



Inserting the limits of s, we get the following bound: 



nn 1 

T> - 

- e TV 

It can be verified that for N » 1, this lower bound is 



T > 



7(1 - 7/4) 



^/N 



(35) 



(36) 



which is a quadratic improvement over the previous bound (32). Also since relation (33) has to be satisfied at every 
instant, the bound (36) is optimal. 

To conclude, in this paper, we have implemented Deutsch's algorithm using adiabatic quantum evolution by a 
Hamiltonian which takes a given initial state to a final state and such that the final state depends on the nature of 
the function /. In particular, if the fimction is constant, the final state is |0) with a very high probability, and if it 
is balanced then the outcome is |1) almost with certainty. Then a measurement on the final state helps to determine 
the nature of /. We have also estimated the required evolution time for a given accuracy of the result. Finally, 
we have generalized the result for the Deutsch-Jozsa problem, using n-qubits, and found that the number of time 
steps required to solve the problem scales as VN, where AT = 2". Although this is a marginal improvement over the 
classically required exponential time (of order N), it docs not match the polynomial time that is achievable using 
standard quantum computational techniques [10]. It would be interesting to compare adiabatic and standard quantum 
computational methods for other algorithms to see whether this difference in computational time is the exception 
rather than the rule. We hope to report on this elsewhere. 
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